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» I ' Abstract. We present a general and novel approach for the reconstruction of 

^ any convex d-dimensional polytope P, assuming knowledge of finitely many of 

its integral moments. In particular, we show that the vertices of an JV-vertex 
convex polytope in K. d can be reconstructed from the knowledge of O(DN) 
axial moments (w.r.t. to an unknown polynomial measure of degree D), in 
d + 1 distinct directions in general position. Our approach is based on the 
collection of moment formulas due to Brion, Lawrence, Khovanskii-Pukhikov, 
and Barvinok that arise in the discrete geometry of polytopes, combined with 
what is variously known as Prony's method, or the Vandermonde factorization 
of finite rank Hankel matrices. 



1. Introduction 

The inverse problem of recognizing an object from its given moments is a funda- 
mental and important problem in both applied and pure mathematics. For example, 
this problem arises quite often in computer tomography, inverse potentials, signal 
processing, and statistics and probability. In computer tomography, for instance, 
the X-ray images of an object can be used to estimate the moments of the underly- 
ing mass distribution, from which one seeks to recover the shape of the object that 
appear on some given images. In gravimetry applications, the measurements of the 
gravitational field can be converted into information concerning the moments, from 
which one seeks to recover the shape of the source of the anomaly. 

The goal of this paper is to present a general and novel approach for the re- 
construction of any convex d-dimensional polytope P, from knowledge of its mo- 
ments. Our approach is quite different from the quadrature-based approach that 
is currently used in the literature. Our starting point is the collection of moment 
formulas, due to Brion-Barvinok-Khovanskii-Lawrcncc-Pukhlikov (in what follows 
referred for brevity as BBaKLP) that arises in the discrete geometry of polytopes, 
and is valid for all dimensions |Bri881 ILaw91i IBar92[ IBar91j . and |BR071 Chap- 
ter 10]. We then set up a matrix equation involving a variable Vandermonde matrix, 
with an associated Hankel matrix whose kernel helps us reconstruct the vertices of 
P. We are also able to reconstruct the vertices of a convex polytope with variable 
density, using similar methods. 
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This new approach permits us to reconstruct exactly the vertices of P, using very 
few moments, relative to the vertex description of P. While the computation of 
integrals over polytopes has received attention recently (see e.g. |BBDL + TT] ). our 
work appears to be the first to provide tools to treat the inverse moment problem 
in general. 

A nice feature of our algorithm is that we do not need to know a priori the 
number of vertices of P, only a rough upper bound for their number. The algorithm 
automatically retrieves the number of vertices of P as the rank of a certain explicit 
Hankel matrix. 

In fact, a surprising corollary is that we only require O(Nd) moments in order 
to reconstruct all of the N vertices of P C M. d . Suppose we are solving the inverse 
moment problem in the context of an unknown density function p. An interesting 
consequence of our algorithm is that even though p is unknown, we may easily 
adapt our algorithm to recover the vertices of P in 0(Nd°d) steps, where d° is an 
upper bound on the degree of p. 

Now suppose we simply solve the direct problem of writing down the moments of 
a given vertex set of a known polytope, with a known polynomial density function 
p. In this direct problem, it would take ( d data to describe the polynomial p 
function, because the space of possible polynomials p has this dimension. However, 
for the inverse problem, where p is unknown, perhaps a counter-intuitive conse- 
quence of our algorithm is that we only require 0(Nd°d) data to recover the vertex 
set of P, which might be smaller than ( d . 

In the existing literature on inverse problems from moments, one immediately 
encounters a sharp distinction between the 2-dimensional case and the general d- 
dimcnsional case, with d > 2. While in the former case a well-known quadrature 
formula allows us to solve the problem exactly for so-called quadrature domains, 
where for the latter case one has to "slice up" the domain of interest into thin 
2-dimensional pieces, solve the resulting 2-dimensional problems, and patch up an 
approximate solution from these 2-dimensional solutions. On the other hand, in 
the recent work of Cuyt et al. [CGMV05 the authors can approximately recover a 
general n-dimcnsional shape by using an interesting property of multi-dimensional 
Padc approximants. 

For z £ M. d and each nonnegative integer j, we define the j-th moment of P with 
respect to the density p by: 

Pj(z) := pj, p (z) := J^(x,zy p(x)dx. 

In this text we restrict ourselves to any density function p which is given by a 
polynomial measure, and which does not vanish on the vertices Vert(P) of P. 

We note that only in the Appendix, when we give proofs of the known BBaKLP 
moment formulas below, we will need to replace the real vector z by a complex 
vector, in order to allow convergence of some Fourier-Laplace transforms of cones, 
but otherwise z will always be a real vector. We say that z is in general position if 
it is chosen at random from the continuous Gaussian distribution on M. d . 

Our main result may be formulated as follows. 
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Main Theorem. Let P C R be a d-dimcnsional polytope with TV vertices, and 
suppose we are only given the data in the form of 0{d p N) moments pj^ p (z), for 
an unknown density p G R[x] of degree d p , and for each of d + 1 vectors z G R d 
in general position. Then the data determines P uniquely, using the following 
algorithm: 

(1) Given 2m — 1 > 2N + 1 moments ci, . . . , c^m—i for z, construct a square 
Hankel matrix H(ci, . . . , C2m-i)- 

(2) Find the vector v = (ao, • ■ • , Om-i, 1, 0, . . . , 0) in Ker(H) with the minimal 
possible M. It turns out that the number of vertices N is in fact equal to 
M. 

(3) The set of roots {xi(z) = (vi,z)|v, G Vcrt(P)} of the polynomial p z (t) = 
t N + X^o 1 Q'iP ^ nen e< l ua l s the set of projections of Vert(P) onto z. 

By contrast with a choice of a general position vector z, we also define, for 
a simple polytope P, a generic vector z G Q d to be a vector that lies in the 
complement of the finite union of hyperplanes which are orthogonal to all of the 
edges of P. For non-simple polytopcs P, we will later extend this definition of a 
generic vector z G Q d , in Section [7] 

We furthermore prove in Section [8] that the vertex set Vert(P) C Q d of any 
rational convex polytope P can be found in polynomial time with a probability 
arbitrary close to 1, from the exact measurements of 0(d p N) moments in carefully 
chosen 2d — 1 random generic directions z G Q d - 

In Section [51 wc indicate how Vert(P) C K d can also be efficiently approximated 
even when the data is noisy. 

One punchline of the proof is that an appropriate scaling of the sequence of the 
moments p-j, P (z) (j = 0, 1, . . . ) for a fixed z is a finite sum of exponential functions, 
and thus satisfies a linear recurrence relation (cf. e.g. |Sta971 Theorem 4.1.1(iii)]). 
Then an application of what is variously known as Prony's method, or Vandermonde 
factorization of a finite rank Hankel matrix (cf. e.g. |BLV97] b allows one to find 
(z,v) for v G Vert(P). As these methods are scattered along quite a number of 
sources, we have chosen to present a self-contained exposition for clarity and for 
ease of efficient implementation. 

Reconstructing Vcrt(P) from the (z,v) is then relatively straightforward, pro- 
vided that wc know these projections for sufficiently many z in general position. 
For the latter, we present an exact procedure as well as a parametric one — the 
latter with the focus being less noise-sensitive. 

The remainder of the paper is organized as follows. In Section [21 we define 
the objects we are dealing with, as well as the appropriate background for ease of 
reading. We also give the known formulas for the moments of simple polytopes. In 
Section [3] we construct a polynomial whose roots correspond to the projections of 
the vertices onto directions z in general positions, by using the moment formulas 
and an associated Hankel matrix. In Section [4] we extend the latter to the case 
of unknown polynomial measures. In Sections [S] and [51 we extend the algorithm 
from Sections [3] and [H which deals with simple polytopes, to all convex polytopcs. 
This completes the proof of the first claim of Main Theorem. In Section [4] we also 
discuss the question of reconstructing p after Vcrt(P) is found. 

In Section [8] we use univariate polynomials to paste together the projections 
retrieved from Section [3] to build up all of the coordinates of each vertex, not just 
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their projections, completing the proof of the second claim of Main Theorem. In 
Appendix lAl we outline some proofs of the known BBaKLP moment formulas from 
Section [21 using Fourier techniques. 

2. Definitions and moment formulas for convex, rational polytopes 

Here we describe an explicit set of formulas for the moments of any convex 
polytope P C R d . We begin with some combinatorial-geometric definitions of 
the objects involved. To fix notation, our convex polytope P will always have N 
vertices. We say that P is simple if each vertex v of P is incident with exactly d 
edges of P. We first treat the case of a simple convex polytope and then later, in 
Section [5] we provide an extension to non-simple convex polytopes. 

There is an elegant and useful formulation, originally due to BBaKLP |Law91j . 
for the moments of any simple polytope in R d , in terms of its vertex and edge 
data. Specifically, let the set of all vertices of P be given by Vert(P). For each 
v G Vcrt(P), we consider a fixed set of vectors, parallel to the edges of P that 
are incident with v, and call these edge vectors wi(v),. . . u>d(v). Geometrically, the 
polyhedral cone generated by the non-negative real span of these edges at v is called 
the tangent cone at v, and is written as K v . For each simple tangent cone K v , 
we let | det K v | be the volume of the parallelepiped formed by the d edge vectors 
wi(v),...,W(i(v). Thus, |detAT v | = \ det(wi(v), . . . ,u>d{v))\, the determinant of 
this parallelepiped. 

The following results of BBaKLP [Law91j give the moments of a simple polytope 
P in terms of the local vertex and tangent cone data that we described above. For 
each integer j > 0, we have 

(1) ^ z ) = f^ E <v,3> J ' +rf D v (z), 

where 



(2) £> v (z) 



v6Vert(P) 

I det K\ 



IlLl( W fc( V )> Z ) 



for each z G R d such that the denominators in D v (z) do not vanish. Moreover, we 
also have the following companion identities: 

(3) 0= E < v ' z >^v(z), 

w6Vert(P) 

for each < j < d — 1. Thus, for example, if z = (1,0, ... , 0), the equations (TTJ) 
and J3]) deal with the first coordinate of each of the vertices v G Vert(P). We 
also note that all of these formulas involve only homogeneous, rational functions of 
z = (zi,...,Zd). 

In the more general case of non-simple polytopes, we may triangulate each tan- 
gent cone into simple cones, thereby getting a slightly more general form of (fTj) 
above, namely: 

(4) ^) = §TdT< £ <v,3) J ' +rf £v(z), 

U vGVert(P) 

where each D v (z) is a rational function that now comes from the non-simple tangent 
cone at v. We note that D v (z) is in fact a sum of the relevant rational functions 
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D v (z) that are associated to each simple cone in the triangulation of the non-simple 
tangent cone K v . 

Throughout the paper, we will mainly work in the context of a continuous domain 
for our choices of admissible z vectors. To be precise, we say that a vector z is in 
general position if none of the following conditions is true: 

(a) z is a zero or a pole of D v (z) , for any v £ Vert(P). 

(b) There exist two vertices vi,v 2 £ Vert(P) such that (vi,z) = (v 2 ,z). 

In the penultimate section, we indicate how to implement our algorithm by transi- 
tioning all of our formulas to a rational context, and picking our z vectors to lie in 
a finite, rational <i-dimcnsional cube. 

The goal here is to reconstruct the polytopc P from a given sequence of moments 
{/•tj(z) j — 1, 2, 3, . . .}. That is, we wish to find an explicit algorithm that locates 
the vertices v of the polytopc P, in terms of the moments of P. To emphasize the 
fact that the moment equations above can be put into matrix form, we define our 
scaled vector of moments by: 

/ x / s f dl(-l) d (l + d)U-l) d kl(-l) d \ 

(5) ( Cl ,...,c fe+1 )= (o,...,0,-^-^Mo,- ^-t-iJH,..., T^rh^k-d) , 

so that the vector c = (ci, . . . , Ck+i) has zeros in the first d coordinates, and scaled 
moments in the last k + 1 — d coordinates. Thus, putting the moment identities ([T]) 
and §3§ above into matrix form, we have: 



(6) 



/ 1 1 

(vi,z) (v 2 ,z) 

(vi,z) 2 (v 2 ,z) 2 



\(vi,z) fc (v 2 ,z) fc 



1 \ 

(vjv,z) 
(vjv,z) 2 

(vjv,z)V 




Recalling that we seek to find the vertices of the convex polytope P with these 
given scaled moments Ci, we answer this question completely, giving an efficient 
algorithm to recover the vertices of such an object P. 

Wc will treat each D Vi (z) as a nonzero constant, which we have not yet discov- 
ered, and each (vj, z) as a variable, for a fixed real vector z £ K d . Moreover, we re- 
alize below that, given our algorithm, only finitely many moments /xi(z), . . . , //m(z) 
are needed in order to completely recover the full vertex set Vert(P), a rather useful 
fact for applications. 

Wc recall that our j'th moment of P is defined, for the uniform measure p = 1 , 

by 



(7) 



(j,j(z) = J (x,z} J dx 



We note that /^o( z ) = Vol(P), the volume of P with respect to the usual Lebcsguc 
measure. 

It is very natural to study moments in this form, because they are "basis-free" 
and also appear as moments of inertia in physical applications. It is worth noting 
that there are other types of moments in the literature, and we mention some 
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connections here. For each integer vector m, we define 



(8) 



Mm = 



Mx, 



with the usual convention that x m = J\i=i x T i an ^ l m l = m i + * • ' + m d- The 
usual application of the binomial theorem gives us a trivial relation between these 
moments: 

k 

mi, . . . ,m d/ 



(9) 



(z,x) = 



E 



m 1 ,...,md- 

roH hmj=t 



In fact, given Vert(P) and n\ m \ (z) as a function of z, we can also compute /x m , as 
following Lemma shows. Its proof is trivial , but it nevertheless offers an interesting 
relation between the moments (J7J and ©• 

Lemma 1. Let m G Z^_, Vert(P), and n\ m \(z) be given. Then 

9l> 



|m|!/u r 



<9z' 



rM|m|(z)- 



3. The inverse moment problem for polytopes - computing 
projections of vert(p) 

In this section we show how, for a given general position vector z, to retrieve 
the projections (v, z) for each vertex v of P, using a certain Hankcl matrix that we 
define below. For the sake of convenience we let Xi = (v;,z), for each 1 < i < N. 
Thus, our goal for this section is to find all Xi, given a number of moments. Due 
to our choice of z, we may assume that xi ^ Xj for i ^ j. From we have 



ci 



\ c k+l. 





/ 1 


1 . 


■ M 




X\ 


x 2 ■ 


■ XN 


(10) 


x\ 


x 2 


T 2 




u 


T . A' 
x 2 


X N/ 



where c is defined by §5§ above 




To streamline notation further, we define a (k 



1) x N Vandermonde matrix ~Vk(x\, . . . , xn), with ij'th entry equal to x* : 









(1 


1 . 


■ M 








X\ 


x 2 ■ 


■ XN 


(11) 


V fc (si,. 


.,X N ) = 


2 

xf 


x 2 


x N 








w 


x 2 


x%J 



We also define a column vector D(z) = (D Vl (z), . . . , D VN (z)) , so that © reads 

Y k {xi 7 . . . ,x N ) ■ D(z) = c. 

We may multiply both sides of (flU)) on the left by a row vector a = (ao , Oi, . . . , afe). 
First, we see that 

A; 

a ■ Vk(xi,x 2 , ■ ■ -,x N ) = {q & {x 1 ),q a (x 2 ), ■ ■ . , fe^Ar)), where q a {t) = ^atf*. 
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Therefore, taking a to be the coefficient vector of the polynomial 

N N-l 

(12) Pz (t) = l[(t-x i )= [] (t-(v,z» 

i=l veVert(P) i=0 

and multiplying (flU)) by a, wc obtain the identity = a • c. Moreover, for each 
< I < k — N wc substitute for a the vector a^ corresponding to t e p z (t), to obtain 
zero in (|10[) . when multiplying on the left by &£. Wc thus obtain k — N + 1 equations 
of the form 

(13) & e ■ c = 0. 

As I increases, the coefficient vector of t e p z (t) gets shifted to the right, and it is 
convenient to capture all of its shifts simultaneously by the to x to Hankel matrix 
H := H(ci, . . . , C2m-i), where we fix m > N + 1, defined by: 



(14) 



H(ci, 



, C 2 r 



Cl 
C2 



C2 
C3 



Cm \ 
Cm+1 

C2OT-1/ 



Theorem 1. The Hankel matrix H has rank iV and its kernel is spanned by the 
to — N linearly independent vectors 



(15) 



at = (0,...,0,a o ,...,aAr-i,l, 0,...,0), < f < m - JV - 1, 



i-N-i-e 



Proof. For each I in the range < t < 2m — 2 — N we use (TT^]). with a vector a^ 
of length 2m — 1. Putting all of these equations together in a more compact form, 
we can write a^H = 0, where now the length of a^ is to. 

It remains to show that the vectors a^ generate the full kernel Ker(H) of H. 
Let b' G Kcr(H). Without loss of generality, there exists an L < N such that 
b' = (bo, . . . , 1,0,..., 0), because we may use the various vectors &£, which lie 
in the kernel of H to get the appropriate zeros in this b' vector. 

Now let us define a number of row vectors of the size 2m — 1: 



(16) b, = (0,...,0,6 ,---,6l-i,1, 0,...,0 ), < £ < 2m - L - 2. 

I times 2ra-2-L-f 

By definition of the Hankel matrix and since to > L + 1, we have b^ • c = 0. 
Consider the polynomial pb(t) = bo + bit + . . . + bi^it 1 "^ 1 + t L corresponding to 
b . 

Taking k = 2m— 2 in (fT0|) . we multiply both sides of (fT0|) on the left by b^. Hence, 
we get hi ■ V2 m -2(si, ■ • ■ , xn) ■ D = 0. Therefore, for every < I < 2m — L — 2 we 
have 

(17) [x[p h (xi), x l N Pb{x N )) • D = 0. 

Combining the first ./V of the latter equations into a matrix form (note that 
N - 1 < 2to - L - 2) we get: 
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/ Pb(xi) 
xlp h (xi) 



\xi 1 pb{xi) 
which can be rewritten as 



pUxn) \ 

XNPb{xN) 

x 2 N Pb{x N ) 



x% 1 pb(x N )J 




( 1 • 

Xl 

xj . 


• 1 \ 

XN 

x N 


(Pb{x x ) 





Pb{x 2 ) ■ 


.. \ 



W- 1 • 


x N j 


V o 





■ • Pb(x N )J 



(D Vl (z)\ 


- 








3 



Since Vjv-i(:Ei, • ■ • , xjv) is invertible, we get 

/ pb(xi)D Vl (z) 



>'Ct 






3 


1- 





\Ph(xN)D VN (z)j 

As L < N and Pb(t) ^ 0, we deduce that x\, . . . , Xn cannot all be roots of Pb{t). 
It remains to mention that D Vi (z) ^ for every 1 < i < N, by the choice of the 
vector z in general position. We therefore arrive at a contradiction. □ 

Once we construct the kernel of H, we will pick a vector (do, . . . , djv-i, 1, 0, . . . , 0) 
in Ker(H), and then define the polynomial p % {t) = a + ait + . . . + aN-it 1 "^ 1 + t N . 
We note that this is the unique vector with the largest number of zeros on the 
right (in algebraic terms, all the remaining vectors in the kernel can be obtained as 
coefficients of polynomials in the principal ideal (p z ) C C[t]). By Theorem [TJ the 
roots Xi (= (vj,z)) of this polynomial are precisely the projections (vj,z) that we 
are seeking. 

(18) p z (t) =a Q + a 1 t + ... + a N - 1 t N - 1 +t N = (t-(v,z». 

veVert(P) 

In summary we have proved the following: 

Theorem 2. Given the moments (|7|) for a direction z G K d in general position, all 
the projections (v, z), v G Vert(P) are the real roots of the univariate polynomial 
p z defined in (fl%|) . 

Finding the kernel of H and then computing the coefficients of p z (t) can be done 
efficiently in polynomial time. After having computed the projections onto z of all 
the vertices, the next step is to find the projections on each of the d coordinates 
of all N vertices of P. However, there is still an inherent ambiguity in this process 
because we will not know from which vertex a specific projection came from. We 
resolve this problem in Section [5] and also in alternative way in Section [5] by using 
univariate representations. 
Remark 3.1. 

(a) An analogue of BBaKLP formula for d = 2 was known for quite a long time 
(see e.g. P. Davis |Dav64j ). and the system of equations corresponding 
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to (fTOj) was solved by what is known as Prony's method, see e.g. Elad, 
Milanfar, and Golub |GMV00| IEMG04] . The solution method described 
above can also be considered as a variation of the Prony's method, 
(b) Importantly, the quantities D Vi (z) play no role for computing the projec- 
tions (vj, z})! This is why we will be able to extend the present methodology 
to general convex polytopcs P (i.e., non necessarily simple). 



4. Polynomial density 

In this section we address the case of non-uniform measures, 
mcnts are now defined as 



That is, our mo- 



(19) 



(x, z) J p(x)dx, 



where the density function p is a homogeneous polynomial of fixed known degree 
d° . We note that, intuitively, if p is not a homogeneous polynomial the change of 
a physical scale (e.g. meters to centimeters) will cause complicated changes in the 
formulas for moments. Therefore, the case of a homogeneous polynomial measure 
is a very natural one, and we begin with this case in order to develop the proper 
formulas for it. We then notice, in the next subsection, that the results for the 
general case of a polytopc with any polynomial density follows exactly the same 
analysis as the case of the homogeneous density. 

To set notation, we let P be a convex polytope with a density function p(x). We 
separate p into its homogeneous polynomial pieces, by writing p(x) = J"^_ Ps(x), 
where p s (x) is a homogeneous polynomial of degree s. We will require the physically 
natural assumption that p(x) > for each x G P, and in fact we will only need the 
assumption p(v) / for v e Vert(P). 

We define Vfc = Vfc((vi,z), . . . , (v^r,z)) = Yk{x\, . . . ,xn), the standard Van- 
dermonde matrix. We further define the Z'th derivative of the Vandermonde matrix, 



namely Vj^ , whose r/'th entry is equal to (i — 1) • (i — 2) 



(t - l)x\ 



(20) 



vP(*i, 



,xn) 




/! 



Jk-i 



(fc-;)! 




/! 



.k-l 



(fc-;)! 



,fc-i 



As mentioned above, we first assume here that p(x) is a homogeneous polynomial 
of degree d°. However, in the following subsection we will discuss how the follow- 
ing formulas also work in the more general case of variable but non-homogeneous 
polynomial density measures. We recall the moment formulas for variable density, 
for a simple polytopc P, from Theorem [S] in the Appendix: 



(21) ^(z) 



j!(-l) d 
(j + d + d°)\ 



y> I _d_ _d_ 



veVert(P) 



dz d 



(v,z) 



j+d+d c 
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where 
(22) 



D v (z) 



dct K v 



n fc =i( w fe( v )' z ) 

and the identity is valid for each z G C d such that the denominators in D v (z) do 
not vanish. In addition, we also have the following companion identities: 

(23) o = £ (v,->^W, 

V 7 veVert(P) 

for each < j < d+d° - 1. 

We now repeat the same procedure of putting the new moment formulas above 
into matrix form, as in ([5]). Here, the definition of the vector c is only slightly 
different, namely: 



(24) 



(Cl, 



,c fe+1 ) = (-l) c 



o o (d + d°y (i + d + d°y 

U ' • ■ • ) U ' 77, MCb Tj Mil • ■ 



fc! • /i fc _ 



d-d" 



0! 



1! 



' (k-d-d°)\ 



We arrive at the following interesting matrix ODE for moments with homoge- 
neous polynomial density: 



(25) 



d 



d 

dz d 



( 1 

(vi,z) 



1 

<V2,Z> 



(vi,z)- (v 2 ,z)' 



\(vi,z) fc (v 2 ,z)' 



1 \ 

(vjv,z) 

(VAT,Z) 2 



(v N7 z) k J 




Cl 



v c fc+l, 



where the differentiation is taken separately for each entry of the vector on the left 
hand side. 

One may check that a single partial derivative of a matrix product obeys the same 
rule as the derivative of a product of two functions, that is -j^{Mi{x) ■ M 2 {x)) = 
£M 1 (x) ■ M 2 {x)) + M x (x) ■ £M 2 (x) 



We compute a partial derivative of Vfe((vx,z), . . . , (vjv,z)): 



dz l 



-V* 



2v«( Vl ,z) 



Vfcv«( Vl ,z) 



fc-1 



V 



(i) 



/v« 



V 



2v«{ Vj v,z) 

fcv 7V ) ( V A r ' Z ) fe_ V 

o \ 





'N 



By repeating the partial derivative in each variable Zi, we arrive at: 
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d_ 



(26) Phrr.---.-5r: 



d_ 

' dz d 



v fc = v fc 



(d°) 



/P(vi) 
p(v 2 ) 



V o o ... ^(vjv); 

Now expanding the matrix ODE formula (|25p . and using the product rule for 
differentiation of matrices, we may write it in the following form: 



d° 






( C1 \ 










j=0 









(27) 



where each entry /| (z) is a rational function of z, and the highest vector term, 
comprised of the rational functions / - d ^ (z) , has the nice form 



(28) 



{f {d ° ] ( Z )\ ( piyi) ° 
> h [Z) \ p(v a ) 



) /^(z) 1 

\A,*(*)y 



V ... ^(vjv)/ 

For the proof of the following theorem we construct a certain vector as follows. 
Define the polynomial 

(d° + l)N-l 

(29) Pz (t)= n (*-<v,z)r +i =^™+ j2 

veVert(P) i=0 

We define the vector to be the coefficient vector of the polynomial t e p z (t) . 

Theorem 3. The Hankcl mxm matrix H, with m > (d° + 1)N + 1 corresponding 
to the moment formulas with variable density, has rank (d° + l)iV, and its kernel 
is spanned by the linearly independent vectors ag. 

Proof. We repeat the procedure that we used in Section [31 using a corresponding 
mxm Hankel matrix and its kernel, but this time the dimension is m > (d°+l)N+l. 
Similarly to the case of uniform density p{x) = - 1, we may again multiply both sizes 
of ([27)) on the left by a row vector ao = (do, a%, . . . , a^-). First, we see that 

a • ...,x N ) = (pi t) (x 1 ) :P ^(x 2 ) 7 . . . ,p®(x N )), 

where Pz \t) is i'th derivative of p z (t). 

Now, multiplying (|27[) by ao, we obtain the identity = ao • c. Similarly, for 
each < £ < k — N(d° + 1), we substitute for ao the vector a^ corresponding to 
t e Pz(t), to obtain = ■ c. Hence the vector ao lies in the kernel of H. 

On the other hand, we now claim that no other vector b different from those 
spanned by a^ could be in Kcr(H). If, contrary to hypothesis, we could find such a 
vector b, we may assume without loss of generality that b = (bo, ■ ■ ■ , h, 1, 0, . . . , 0), 
with I < (d° + l)iV — 1, by reducing it with appropriate linear combinations of a.£. 

Recall that c = (a,..., c^+i), where k > 2m — 2 > 2(d° + l)N. Let us consider 
polynomial Ph(t) with coefficients of b. Now let hi corresponds to the polynomials 
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t e Pb{t), for < I < m. Then we have • c = 0, because b is in the Ker(H) and 
each vector hi has the same entries as b only shifted by t to the right. 

Since a degree of Pb(t) is smaller than that of p z (t), we have p z (t) \ Ph(t)- 
Therefore, there exists v £ Vert(P) such that (t - (z,v)) d ° +1 \ Pb(t). Without 
loss of generality, we may assume that v = Vi. We now construct a polynomial 
q(t) by multiplying Pb(t) by sufficiently many linear factors of the form (t — (z, v)), 
where v varies over all of the vertices of Vert(P). We will treat the particular 
vertex vi differently, by multiplying by a slightly different power of (t — (z, Vi)), 
to insure that a certain derivative, explicated below, does not vanish at xi, thus 
giving us a nonzero vector in the kernel of H. The desired polynomial q(t) satisfies 
the following properties: 

(1) Pb(t)\q(t). 

(2) {t-{*,v x )) d °^\q{t). 

(3) (t-fcvx))* \q(t). 

(4) (t - (z,v)) d ° +1 | q(t), for Vv G Vert(P) : v ^ v x . 

(5) deg(g) < dcg(p) + N(d° + 1). 

We now write the coefficients of polynomial q{t) as a vector b°. Next, we multiply 
((27)) on each side by the row vector b°. 



(30) 



E b °- V 



W 

k 













Vk+ij 



The vector b° may be represented as a linear combination of vectors hi, where 
< I < m. Therefore, we get b° • c = 0. 

On the other hand, since ~ x i) d we have b° • V| = for each 

< i < d°. Then 



h°.vP = (q( d °\ Xl ),..., q ^(x N j), 

where xj = (z,Vj). We have q( d °\xi) = 7 =^ 0, by property ([2]), and q( d \xi) = 
for each 2 < i < N ', because n^2(^ — x i) d " +1 \<l{t)- Therefore, we get 



(31) 



(7,0,. ..,0) 



/>(vi) 
p(v 2 ) 



V 







\ 





p(vat)/ 



'D Vl (zY 



= 0. 



Thus 7 • p(vi) • £> Vl (z) = h° ■ c = 0, where none of the quantities 7, /o(vi) and 
D Vl (z) is zero, so that we have arrived at a contradiction. □ 

Therefore, we have proved the following result, the analogue of Theorem [2] for 
the homogeneous polynomial density case. 

Theorem 4. Given moments ([T9| for a direction z 6 K d in general position and 
where p is a unknown homogeneous polynomial of degree d°, all projections (v,z), 
v G Vert(P), are the real roots of the univariate polynomial p z defined in 
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4.1. Non-homogeneous measure. We start with the moment formulas for a 
polytope with variable, but homogeneous density, namely (|54]l. Now we let d° 
be the maximal degree of the monomials of p(x). Then the formula (|21[) can be 
rewritten as follows. 



d° oo , v j oo 

(32) x>.(v.)t- e E^c-^v^-^Ef^'- 

s=0 vGVcrt(P) j=0 J ' j=0 7 ' 

Following the same reasoning that was used for the homogeneous variable density 
case, we first collect all the coefficients of p~ d ~ d on both sides of (pl2"T) . to get: 



(33)$>(V.) £ " 1} ■ • • d ° + S + 1) (v, zr*°+« -^(z) = 3±1, 

s=0 vGVert(P) ^ ' ^ ' 

where c J+ i = (— l) d J ^ J d l d do d ^ , as in the formula (f2"4"|). Next, we put everything into 
a matrix form and get 



(34) 



5=0 



p s (Vz) 



V 



(xi, . . . ,x d ) ■ D 



\CfeH 



The latter matrix ODE can be brought into the same form as (|2"T)l . with exactly 
the same coefficient of (xi, . . . , xn) that appears in (|2"5|) . Therefore, our method 
works for general polynomial density measures as well, with precisely the same 
algorithm. 



5. General convex polytopes 

In the previous discussion we considered only simple polytopes, because the 
BBaKLP formula takes a particularly nice simple form when P is a simple polytope. 
However, it is natural to extend our approach to non-simple polytopes. Indeed, 
it is always possible to triangulate P, that is decompose P into a union of non 
overlapping simplices, without adding any extra vertices (See, for example, |BR071 
Theorem 3.1]). 

We now fix one such triangulation of P, and denote it by T(P). We may then 
rewrite the formula for each moment /Ltj(z) as follows. 



(35) ^(z)= / (x,z)^x= E / (x,z)^x. 

Triangulating the general convex polytope P into simplices, we reduce the general 
moment problem to the moment problem for each simplex A of the triangulation. 
Although triangulations may be expensive to construct in practice, we only need to 
consider a theoretical non- vanishing result, given in Lemma [5] below, for any such 
triangulation. Given such a triangulation, we may then apply the formulas ([1} and 
§5$ to each of the simplices A in the equation above: 
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c » = E E Wa,(a, Z ) 

AeT(P) vSVcrt(A) 

= E < v > z >' E ^v(A, Z ), 

vGVert(P) Ae{T(P) | vGVcrt(A)} 

where we have interchanged the order of summation in the last equality above. We 
now define D v (z) for this fixed triangulation T(P) by: 

(36) D v (z):= E D v (A,z) 

Ae{T(P)|veVcrt(A)} 

Then we have 

(37) Cj (z)= ( v - z )^v(z). 

veVert(P) 

This gives us 

(38) {j + d f- 1) \ j (z)= E (v ; z)^ + ^ v (z). 

vGVert(P) 

Note, that in Section[3]we never used the explicit formula for D v (z). The only fact 
we exploited was that D v (z) ^ for a general position vector z. Therefore, we 
can apply the same approach for non-simple polytopes, if we are able to prove that 
D v (z) 7^ for a general position vector z. 

Lemma 2. For any vertex v G Vert(P), any fixed triangulation T(P) and a general 
position vector z we have D v (z) / 0. 

Proof. We begin by noting that D v (z) is a finite linear combination of rational 
functions of z. In fact, according to the Lemma 8.3 and Chapter 9 of |Bar08j . 
D v (z) is a rational function that is the analytic continuation, in z, of the function 



i Xv -v(z) = / e < z '*- v W, 

J K v 

when this integral converges. We define the dual cone to K v — v as follows: K* := 
{y G M. d (y, x) < 0, for all x G K v — v}. Indeed, the latter integral converges 
for all z lying in the interior of the dual cone K*. Since K v is a tangent cone of a 
convex polytope, the dual cone K* is non-empty. Clearly e^ z ^ is positive for all 
x G A' v — v, if z G K*. We obtain the result that l^ v _ v (z) > for all such z, 
and we may therefore conclude that the analytic continuation of l^- v _ v (z) cannot 
vanish. □ 
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6. An exact algorithm 

In Section [3] we have learned how to find the projections of vertices of P onto 
a general position axis z. A short summary of the procedure for such a randomly 
picked z € R d is as follows: 

(1) Given 2m — 1 > 2N + 1 moments c±, . . . , C2 m — l f° r z > construct 
a square Hankel matrix H(ci, . . . , c% m — l)- 

(2) Find the vector u = (ao, . . . , a,M-i, 1, 0, . . . , 0) in Ker(H) 
with the minimal possible M. It turns out that the number of 
vertices N = M. 

(3) The set of roots {xi(z) = (vj,z)|vj £ Vert(P)} of polynomial 
Pz{t) = ao + ait + . . . + aiy-it 1 ^^ 1 + t N then equals the set of 
projections of Vert(P) onto z. 

Algorithm 1: Computing projections. 

Remark 6.1. Note that N is an essential part of the input. One cannot rule out 
existence of another polytope P' with Vert(P') > N and the same moments, up 
to certain degree. 

Remark 6.2. If we work in the context of exact measurements, with rational 
vertices and rational choices of z vectors, then p z has only rational roots. In this 
rational context, we may analyze the complexity issues involved by using the LLL- 
algorithm due to Lenstra, Lenstra, and Lovasz |LLL82| . because now the rational 
roots of p z can be found in time which is polynomial in N and in the bitsize of 
Vert(P). 

In this section, we describe below an exact algorithm to compute Vert(P) that 
runs in polynomial time given the latter assumptions. When the roots of p z are 
not available exactly, the algorithm still works, producing approximate results. 
However, it seems nontrivial to control the precision of root-finding, as we need 
to find the roots of d univariate polynomials. In Section [5] we present a different 
procedure, where, in contrast, roots of only one polynomial parametrize Vcrt(P), 
and which conceivably is more robust against numerical errors. 

We use the assumption that z is in general position (it suffices to require that z 
is not perpendicular to the lines uv, for u,v £ Vert(P)) to maintain bijectivity of 
projection onto z, as well as to avoid division by zero in the terms D Vi (z). Choosing 
z at random from the Guassian distribution on K d , we get a z in general position 
with probability 1. Further, to reconstruct the locations of Vert(P) given the pro- 
jections of vertices on a number of axes we match all projections of the same vertex 
as follows. 

• Take d linearly independent vectors zi, . . . each chosen in general po- 
sition. 

• For every 2 < i < d match projections of Vert(P) onto z, with projections 
onto Zi. 

(1) Pick a general position vector z = azi + f3zi in the plane generated by 
zi and z^ 

(2) Compute the coefficients of the polynomial p z (t) using extra 2N + 1 
moments in direction z. 
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(3) For each pair of projections Xj(zi), Xfe(zj) onto zi and z, match them 
whenever p z (axj + (3xk) = 0, for 1 < j, k < N. 

(4) With probability 1 all vertices will be matched correctly, that is Xk(zi) 
is matched with Xk(zi)- 

• For each 1 < k < N reconstruct v k 6 Vert(P) from its projections Xk(zi) 
for 1 < i < d. 

Indeed, the degree N polynomial 

P*(t) = H (t - ax k (zi) - /3x k (zi)) 

k 

has N distinct roots. We evaluate it at the N 2 values QXj(zi) + j3xi(zi). With 
probability 1, by choice of a and /3, p z will only vanish when Xj(zi) and xg(zi) 
correspond to the pro jections of the same vertex. 

(In fact, this part is easy to de-randomize: fixing a = 1 and choosing more that 
N 3 different values of /? gives one a "good" pair a, f3.) 

Note that in total we have used (2c? — 1)(2A + 1 — d) distinct moments, while 
the description of vertices of P requires d ■ N real numbers. That is, our procedure 
is quite frugal in terms of the moment's measurements. 

As claimed in Main Theorem, we can still improve on the latter (albeit the 
corresponding procedure is not polynomial time any more). Indeed, we only have 
moments for d + 1 directions Zi,. . . , z c i, z = . otjZj in general position, we can 

still carry out a similar procedure, although one would need to compute test 
values (for all the possible d-fold matchings) of 

d 

k j=i 

7. AN ANALYSIS OF OUR ALGORITHM IN THE RATIONAL CASE 

In Section [5]we described our algorithm under the global assumption that each 
direction z is chosen at random from the continuous domain K . , thus getting a gen- 
eral position vector z with probability 1. However, in any practical implementation, 
all the coordinates of z have to be rational numbers with bounded denominators 
and numerators. In this case the probability that z does not lie in general position 
will be strictly smaller than 1. In what follows we describe a way to pick our z- 
directions and argue that the probability for choosing a "bad set" of z-directions 
(which are not in general position) is indeed small. 

We will always pick our z vectors to be rational vectors, with denominator equal 
to r, and lying in the unit cube [0, l] d . If we knew the vertex description of a simple 
polytope P, we would only need to make sure that z lies in the complement of the 
finite union of hyperplancs that arc orthogonal to all lines between any two vertices 
of P. We call such a rational z a generic vector. The probability of picking such a 
generic z tends to 1 as r — > oo. 

We now extend the definition of a generic vector z to a non-simple polytope P. 
In this case, in addition to our previous restriction that z is not orthogonal to any 
line between vertices of P, in particular to the edges of P, it might occur that z 
is a zero of the rational function D v (z), defined by (|36| in Section O and we need 
to avoid such a choice of z. Hence we define a generic vector in the general case of 
non-simple polytopes to be a vector that is simultaneously not orthogonal to any 



INVERSE MOMENT PROBLEM FOR CONVEX POLYTOPES 



IT 



line between vertices of P, and also not a zero of any rational function D v (z). In 
particular, we shall avoid zeros and poles of the complex function D v (z) in z. 

In what follows, we refer to the algorithm of Section [5] By the Schwartz-Zippel 
Lemma |Sch801 |Zip79| IDL78] , we have an upper bound for the probability that the 
numerator and denominator of the multivariable rational function D v (z) vanishes 
for a random rational z E [0,l] d , where z has denominator r. In fact, by our 
construction, we have r d such rational vectors z, and the Schwartz-Zippel Lemma 
tells us the following: for sufficiently large prime r Prob[z is a zero of D v (z)] < ^ 
and similarly Prob[z is a pole of D v (z)] < ^. Indeed, both the numerator and 
the denominator of D v (z) are homogeneous polynomials in d variables Z\ . . . , Zj, of 
degree at most N with integer coefficients; none of these polynomials vanish when 
taken over the finite field ¥ r , for all sufficiently large r. 

We remark that our algorithm picks either arbitrary generic vectors (we pick 
them uniformly at random from the rational unit cube), or takes an integer linear 
combination of two independent random vectors. In the former case by taking r of 
order 2 polv ( N,d ) one can make the above probabilities for all D v (z) to be negligibly 
small. In the latter case, we need to be more careful, as the sum of two random 
vectors uniformly distributed over the rational unit cube is no longer a random 
vector distributed uniformly over the unit cube. However, we may now consider 
the vector az\ + /3zi, as well as the numerator and denominator of D v (z), over the 
finite field F r . We note that, once we fix < a < r and < (3 < r, the linear 
combination of two independent, uniformly distributed vectors, namely az\ + /3zj, 
is again uniformly distributed over Ff. 

Therefore, we may assume that each particular direction z that appeared in the 
algorithm [1] is generic with a very high probability. On the other hand, a generic 
vector z = az\ + /3zj in the plane spanned by zi,Zj, matches the set of projections 
onto zi and the set of projections onto Zj uniquely at very high probability. Indeed, 
given the projection onto Zi and Zj there are ./V 2 possible projections of Vert(P) 
onto the plane spanned by Zi and Zj and at most ./V 4 different lines between these 
points. In other words, there are altogether at most iV 4 directions that do not help 
us match projections onto zi and z^. In the algorithm we pick one of the r distinct 
directions for z = azi + /?Zj for any fixed a. Thus the chance that our algorithm 
did make a mistake in a particular step is negligibly small. 

8. Univariate representations for Vert(P) 

In this section, we present an alternative procedure, that is conceivably more 
robust than the algorithm in Section|Hl where given a finite collection of projections 
of the vertices, we presented an exact procedure to reconstruct them. That is, we 
were given some data described in Algorithm [T] assuming that Vert(P) C Q and 
the measurements are exact. When at least one of the latter assumptions does 
not hold, the polynomial p z , whose roots are projections of Vert(P), may not have 
rational roots. Even its coefficients might be known only approximately Thus it 
might be hard to control numerical errors. 

We construct univariate representations (see e.g. [BPR 03 ) of v G Vert(P). The 
latter are typically used to compute solutions of systems of multivariate polynomial 
equations — here this appears to be the first use of these representations for purposes 
other than solving systems of polynomial equations. That is, we will express the 
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coordinates of v £ Vert(P) as univariate rational functions of where •& is a root 
of p a (t) in (fT8|. 

Wc introduce bivariate polynomials / a b £ K[s,i] denned by: 

(39) (s,t) H./ ab (*,t) = J] (t-(v,a+«b», a,bGK d . 

vGVert(P) 

Upon transitioning to rational vectors a and b, generic in the sense of Section [3 
and with a ^ b, we can compute the coefficients of / a t>(s,i) by interpolating, 
with respect to s, the coefficients of the polynomials / a b(s,i) = p a +sb(£)j with 
s = 0, 1, . . . , N, and p a + s h in (fl"8]) computed using Theorem [TJ Define 

(40) .g a b(t) := ^ | a =o ■ 

Then 

ff.b(*)=- e < v ' b ) n 

vGVcrt(P) v^uGVert(P) 

In particular for w £ Vert(P) one obtains 

,9ab((w,a)) = - ^ (v,b) (w-u,a) 

vSVert(P) v^ueVort(P) 

= -(w,b) (w-u,a). 

w^uGVort(P) 

On the other hand, for p a in (|18p . its derivative p' a reads 

p*®= e n (*-<«>*» 

vGVort(P) v^uGVert(P) 

and thus 

p4(( w - a ))= E II (w-u,a)= Yl (w-u,a). 

vGVcrt(P) v#u£Vcrt(P) w^ueVert(P) 

Hence 

, ,v 3ab«W,a}) 3ab(#) , q , , q s n 

w ' b =^77 rv = , , a x , f or some i? s.t. p a (tf = 0. 

In particular, assuming that a set of basis vectors ex, . . . , e^ of K d are generic, we 
obtain 

Theorem 5. The set of vertices of P is given by 

(41) Vert ^) = {(^>--->^)l each * B.t. = o} , 

provided that a, ex,...,e<j € M. d are 'sufficiently general' w.r.t. P - that is, the 
polynomial p a (t) from (fT5|) and the polynomials g ae (t) from (|40[) have no multiple 
root. □ 

We remark that the assumption of being "sufficiently general" in Theorem [5] 
is equivalent to the fact that each of the vectors a, ex , • • • , e^ does not lie in the 
discriminant varieties of the polynomial p a {t) and the set of polynomials g aej (t). 

Assuming that computation is done with arbitrary precision, the vertices of P 
can be obtained by evaluating the vectors of rational functions in i) at the roots 
of p a , as in (|4"Tj) . Therefore, we have transformed the delicate computations of the 
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roots of the polynomials p z for all projections onto a number of axis vectors z, into 
just one calculation given by (|4Tj). 

We note that here we need to use 0(dN 2 ) moments, which is typically much less 
frugal than the method of Section [6] which only uses O(dN) of them. 

Remark 8.1. A similar computation of the univariate representation can be carried 
out even without the genericity assumptions, when the corresponding univariate 
polynomials have multiple roots. See |GP05j for details. 

9. An application to physics 

Here we discuss an application of our results to a classical problem of mathe- 
matical physics — reconstruction of an object from the potential of a field that it 
creates. For concreteness, we limit ourselves to the 3-dimensional potential of the 
gravitational field. The potential function u(x) := u(x\,X2, £3) of the gravitational 
field F(x) is defined by 

F{x) = Vu{x). 

In turn, for a body TcK 3 with density p(x) the potential is given by 

u(x) = ( .. ^ ^ .. dt, for any x T. 
Jt \\ x ~ t\\ 

A typical physics problem is to reconstruct T and p from u, i.e. from the measure- 
ments of u. That is, we can assume that ||x — i|| = Y] fa(x)t a is an expansion 
in a Taylor series w.r.t. t = {ti,t2,ta), and the f a (x) depend upon x only. Then 
the expansion 

u{x) = V f a (x) / t a p(t)dt, for any x <£ T 

encodes information of the moments J T t a p(t)dt of the measure p(t) supported 
on T. Thus reconstructing T and p from u is an inverse moment problem. For 
instance, when p is a polynomial and T is a polytope, the approach described in 
this paper can be applied to this inverse potential problem and will provide an 
exact reconstruction. 

Acknowledgments. We thank the referee for very useful suggestions, which in- 
deed improved the text. 
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Appendix A. Proof of the BBaKLP identities for moments of 

POLYTOPES 

Here we recall the proof of the moment formulas of BBaKLP, as well as some use- 
ful facts that arise in combinatorial geometry, and which we have used in the present 
paper, but which may not be well-known yet to the mathematical community at 
large. Although some of the proofs here may not be completely self-contained, they 
give the reader the proper background for understanding where the moment for- 
mulas come from, and the tools that are used for handling them. For more detailed 
proofs of some of these results, the reader may consult the book jBar08j . as well as 
the book [BR071 Corollary 11.9]. We begin with a very useful geometric identity, 
which has an inclusion-exclusion structure, due to Brianchon and Gram. 

We let lp(x) denote the indicator function of any convex polytopc P. For any 
d-dimcnsional convex polytope P, we have the following Brianchon-Gram identity: 

(42) 1p(x)= £(-l) dim(F) l* F M, 

PCP 

valid for all x G R d . Here we are using the tangent cone Kp at each face F of P. 
Lemma 3. 

ip(x) = y ^W. 

ueVert(P) 

Proof. We simply take the Fourier-Laplace transform of both sides of the Brianchon- 
Gram identity above, and we recall that it is defined by /(z) := / Rd /(x)e^ x,z ^dx, 
valid for all z £ C d for which the integral converges. By definition, we have 



ip(z) = J e< x ' z >dx, 



p 

the Fourier-Laplace transform of the indicator function of P. It turns out that we 
may define the Fourier-Laplace transform 1k f (x) = 0, for any tangent cone Kp 
which contains a line (isomorphic to K 1 ). Since all tangent cones Kp contain a 
line, except for the vertex tangent cones, we are left only with the Fourier-Laplace 
transforms of the vertex tangent cones. Precisely, we get: 

i P (x)= Y W*)- □ 

u6Vert(P) 

Using the theory of valuations, one can make the proof of the former Lemma 
more rigorous (see |Bar08j ). However, for the purposes of this appendix, it is not 
necessary to consider the subtle issues of convergence that arise here. 

Lemma 4. Let K v be a vertex tangent cone of a simple polytope P. Then 

P ( v > z > Het K 



rnU( w fc( v )> z ) 



for all z E C d such that the denominator does not vanish. 

Proof. The main idea here is to use the fact that there is a linear transformation 
that maps the simple tangent cone K v bijectively onto the positive orthant K ort h '■= 
{(xi, . . . , Xd) S M. d Xj > 0}. To be explicit, let K v — v := Kq be the translated 
copy of our tangent cone K v , so that the vertex of K lies at the origin. Let M 
be the invertible matrix whose columns are the d linearly independent edge vectors 
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Wfc(v) of K v . Then the linear transformation T : K or th — > K v — v, defined 
by T(x) = Mx, gives us the desired bijection from the positive orthant onto the 
translated tangent cone K v — v := Kq. Now we use the explicit computation for 
the Fourier-Laplace transform of the positive orthant K ort h, namely: 



d d 



i^w=f[i«>ofe)=H) d n - 



3=1 3=1 

Finally, the standard Fourier identity (/ o T)(z) = | dctT|/(T'z), valid for any 
invcrtible linear transformation T, allows us to finish the computation: 

iif v ( z ) = lKo+v(zs) 

= e^\detM\i Korth (M^) 



e^(-l) d detK v f[( / ) 

M V( w fe( v )> z )/ 



3=1 

□ 

Theorem 6. Let P be a simple convex polytope. An explicit formula for the 
Fourier-Laplace transform of P is given by: 



( 43 ) / e< x < z >dx=(-l) ^ — d 



^ e< v - z >detif v 

rt (P) nLi( w fe( v )' z ) 



for all z that are not orthogonal to any edge of P. 

Proof. From Lemma |3l we know that the Fourier-Laplace transform of P is given 
by the sum of the Fourier-Laplace transforms of the vertex tangent cones K v , over 
all vertices v of P. Using Lemma 2] to rewrite the Fourier-Laplace transform of 
each vertex tangent cone explicitly, we are done. □ 

Theorem 7. For any convex polytope P and any polynomial p G R[x], there exist 
rational functions q v (z) such that 

(44) f e< x ^p(x)dx= e (v ' Z> 9v(z), 

P v6Vert(P) 

for all z such that the function e^ v ' z ^ q v (z) is analytic at z. 

Proof. We may first employ the fact that every convex polytope P has a trian- 
gulation into some M simplices Aj, with no new vertices. We therefore have 
ip(z) = X)i=ilAi( z ), because the d-dimensional Fourier transform vanishes on 
all of the lower-dimensional intersections of the various simplices Aj. We observe 
that 

(45) f.™^-,(jL,...,£.y.*«^ 
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because due to the compactness of P, differentiation under the integral sign is valid. 
Thus 

Now by Theorem [SI applied to each simple polytope Aj, we finally have 

,47) / e<->p(x)* -,( ») t(-V* E f^TTTT - 
Jp i=i veVort(A,) n fe =i(wfc(v),z) 

giving us the desired conclusion upon applying the differential operator to each 
rational function. □ 

We recall from the introduction that the basis-free moments for uniform density 
were defined by 



W( z ) ; = J Jjt,iY dx. 



The following set of moment formulas can also be found in |Bri881 Section 3.2], 
as well as in |BR071 Section 10.3]. 

Theorem 8. (Moments Formula for uniform density) Given a simple polytope P, 
with uniform density p = 1, we have the moment formulas: 

(48) ^.( Z ) = &^ £ <v,z)^ v (z), 
for each integer j > 0, where 

(49) D V ( Z ) := J dCtAv| , 

rL=i( w fc( v )> z ) 

for each z s C d such that the denominators in -D v (z) do not vanish. Moreover, we 
also have the following companion identities: 

(50) 0= ]T (v, Z )^ v ( Z ), 

i>6Vert(P) 

for each < j < d — 1. 

Proof. We begin with the explicit identity for the Fourier-Laplace transform of any 
convex polytope, namely (|43p . and we replace z by tz, where t > is now treated 
as a real variable: 



e*< v ' z > detA'v 



Now we expand both sides in their Laurent series about t = 0, and equate the 
coefficient of V on both sides to obtain the desired moment identities. □ 
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Theorem 9. (Moments Formula for polynomial density and any convex polytope) 
Suppose we have a homogeneous polynomial density function p(x), of degree d°, 
defined over any convex polytope P. For each integer j > 0, we have the density 
moments formulas 

<*> m^^£ 6v5(4/ (a,..,a) (v , z , WDv(z) , 

where 

(52) D v (z) := J dCtXv| . 

n fe =i( w fc( v ) 5 z ) 

These identities are valid for each z € C d such that the denominators in D v (z) do 
not vanish. In addition, we also have the following companion identities: 

(53) = p(A,.. J2 <v,.>^(.), 

V 7 i=l vGVert(Ai) 

for each < j < d+d° - 1. 



Proof. We begin with (j47[) . and replace z by te, for any fixed t > 0. Again, 
expanding both sides in their Laurent expansions about t = gives us: 

J=0 J ' v a/ i=\ vGVort(A,) J=0 J ' 



v 7 i=l veVert(A«) J=0 J 

We now equate the coefficient of for each j > 0, on both sides of the former 
identity (|54[) . to obtain the desired moment formulas for variable density: 

jl(-l) d ^ f 5 5 



Moreover, we also obtain the desired companion identities (|53j) . by equating the 
first d + d° coefficients of (54]), for each < j < d + d° - 1. □ 



